Ramsey fringes in a Bose-Einstein condensate between atoms and molecules 
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effectively to a channel which is energetically closed. In a 
system of two rubidium atoms, the position of this bound 
state depends on magnetic field, due to hyperfine and 
Zeeman interactions. The Feshbach resonance gives rise 
to a dispersive behavior of the scattering length, that can 
be formulated accurately as 
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In a recent experiment, a Feshbach scattering resonance was exploited to observe Ramsey fringes 
in a 85 Rb Bose-Einstein condensate. The oscillation frequency corresponded to the binding energy 
of the molecular state. We show that the observations are remarkably consistent with predictions of 
a resonance field theory in which the fringes arise from oscillations between atoms and molecules. 

PACS numbers: 03.75.Fi,67.60.-g,74.20.-z 

Interest in the physics of ultracold molecules has been 
growing considerably in the past few years One 
of the goals has been the creation of a molecular Bose- 
Einstein condensate (BEC). In a recent experiment at 
JILA, performed by Donley et al. ||, a coherence be- 
tween atoms and molecules was demonstrated in a BEC 
of 85 Rb atoms. In this experiment, two magnetic field 
pulses were applied to the condensate, and oscillations 
observed in the population of two spatially distinguish- 
able components as a function of the pulse separation 
time. One of the components, the 'remnant' atoms, had 
a similar spatial profile to the original BEC. This was 
in contrast to a second component, the 'burst' atoms, 
which moved away from the remnant having gained a 
considerable amount of energy. In addition, there was a 
missing fraction of atoms unaccounted for after the pulse 
sequence. Significantly, the frequency of the population 
oscillations corresponded to the binding energy of the 
highest lying molecular state. 

This observation followed other remarkable experi- 
ments carried out in the same group. The collapse of 
a BEC was studied when a Feshbach resonance was used 
to create a large negative scattering length [Q. Because 
of the rather violent destruction of the collapsing con- 
densate, this effect was dubbed a 'Bosenova' in analogy 
to a supernova explosion. In a precursor experiment to 
the one we consider here, a single strong-coupling field 
pulse was applied, and there remnant and burst atoms 
were also observed j|. 

In this letter, we describe the most recent JILA experi- 
ment by a resonance field theory that has been developed 
over the last two years HQ. This mean-field theory of 
dilute atomic gases goes beyond the level of the Gross- 
Pitaevskii equation ||[l0[ to include the essential pairing 
physics necessary to describe resonances in the two-body 
scattering. In the case of a Fermi gas it was applied to de- 
scribe superfluidity close to a Feshbach resonance. Here, 
however, we are able for the first time to compare this 
approach with experimental data, and we obtain remark- 
able agreement with the observations. 

A critical ingredient of the theory is the description of 
the Feshbach resonance. Such a resonance arises when 
a bound state (molecular state) lies near the threshold 
of the collision continuum. This bound state belongs 
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Here db g is the background scattering length, AB the 
width of the resonance, and Bq the field value where the 
scattering length is infinity. The magnetic field can be 
easily converted into an energy detuning of the molecular 
state from threshold by the relation i/q = (B — Bq)A^l, 
where Afi is the difference in magnetic moments of the 
energetically open and closed channels. In our theory, 
we treat the closed channel explicitly by using molecular 
field operators and ascribing a coupling to the continuum 
for this molecular state. 

The Hamiltonian for the resonance system is given by 

H= Jd 3 x ($l(x)H a (x)# a (x) + ^ m (x)H m (x)vh m (x)) 
1 



d 3 x x d 3 x 2 



^l(xi)ipt(x2)V(x 1 - x 2 )-ip a (x 2 )ip a (x 1 ) 



+ (^(^r^ )ff(a;i - s 2 )&(* a )&(*i) + H.c.)] , (2) 

where the field operators ^ a {x) and ip^x) create an 
atom or molecule at position x, and H.c. denotes the 
Hermitian conjugate. The free Hamiltonians H a (x) = 
-fi 2 V 2 /2m and H m (x) = -% 2 V 2 /im + v for atoms with 
mass m and molecules with mass 2m include the detun- 
ing v. Atom-molecule collisions and molecule-molecule 
collisions give higher order corrections. The potential 
terms V(x\ — x 2 ) and g(x\ — x 2 ) have to be chosen 
such that both the scattering physics and the molecular 
binding energies are correctly described. We verify this 
by noting that the scattering equations are included in 
this resonance mean- field theory M. Although unneces- 
sary here, it is important from a fundamental perspective 
that the above Hamiltonian can be further generalized by 
the inclusion of more scattering resonances to systemati- 
cally improve the description of the two-body scattering 
physics to any desired precision. 
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Starting from this Hamiltonian we obtain the Hartree- 
Fock-Bogoliubov (HFB) equations of motion. We define 
expectation values for the atomic and molecular conden- 
sates (j) a (x) = (-0 Q ) and <j>m{x) — ($ m ). Moreover, we 
define a 2 x 2 density matrix for the fluctuating compo- 
nents of the atomic field operators x a — ip a — (V> n ), that 
describe the noncondensed atoms: 



G(x,y) 



(xUy)xa( x )) (xa(y)x a ( x )) 
(xl(y)xU x )) (xJy)xU x )) 



(3) 



The elements of this matrix can be given in terms of 
the normal density GN( x ,y) = {x\{y)x a i x )) an d the 
anomalous density G A (x,y) = (x a (y)x a ( x )) [§■ 

To begin with, we consider a gas which is homogeneous 
and isotropic, which results in a translationally invari- 
ant system. This implies that the single particle fields 
are constant in space, and that the two-particle fields 
depend on the coordinate difference r = \x — y\ only. 
We substitute local interactions for the potential terms: 
V{r) — VS(r) and g(r) = gS(r), where V and g are con- 
stants. The HFB equations for the density matrix Q and 
for the condensate fields <p a and (j> m are obtained from 
the Heisenberg equations for the field operators, having 
taken expectation values and applied Wick's theorem: 
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dG N {r) 

dt 
dG A {r) 
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2Re [V(4> 2 a + G A (0))G* A (r) + gcf> m G A (r)} , (6) 
V2 -G A {r) + 4V(\^ a \ 2 + G N (0))G A (r) 



to be appropriate zero range limits of nonlocal poten- 
tials (e.g. square well potentials), and the properties of 
these potentials can then be chosen such that the mi- 
croscopic low energy two-body physics around a Fesh- 
bach resonance is correctly described. This renormaliza- 
tion procedure amounts to replacing the coupling con- 
stants in the Hamiltonian by parameters which depend 
on K. The required parameters can be concisely summa- 
rized by the complete relations V = TU, g = Tgo, and 
v = i/Q+aggo/2, where U — 4nh 2 ab E /m, V = (1 — all)^ 1 , 
a = mK/ (2n 2 h 2 ), and go is determined from the field de- 
pendence of the binding energy as we now explain. 

A big advantage for this system is that the rubidium 
two-body interactions are extremely well known. Typi- 
cal scattering lengths can be calculated to at least the 
1% level [[ll]. For this experiment, where the fringe fre- 
quency is determined by the energy of the highest bound 
state, accurate knowledge of the binding energy is cru- 
cial. In Fig. we show the binding energy as a function 
of magnetic field using a full coupled channels calcula- 
tion [EJI . This calculation consists of all relevant rubid- 
ium two-body interactions including hyperfine and Zee- 
man interactions and electron-spin dependent singlet and 
triplet potentials. In order to reproduce exactly the re- 
sult of the full coupled-channels calculation, we would 
need to go beyond the single resonance formulation we 
have presented. However, in spite of the complexity of 
the real rubidium system, the field theory with one reso- 
nance state allows a remarkably good approximation to 
the binding energy over the field range of interest. This 
comparison is shown in Fig. |l| and determines the value 
of go used in our following simulations. 
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with \i the reduced mass. This is the complete closed 
set of equations to be dynamically solved. As empha- 
sized earlier, the binary collision physics encapsulated in 
the HFB equations is extracted by setting the density 
dependent shifts to zero 
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VS(r) 



V(r)+gS(r) 
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where V(r) = <j> 2 + G A (r) is the total pairing field, leav- 
ing the two coupled two-body scattering equations with 
contact interactions ||. 

The local potentials in Eqns. give rise to an 

ultraviolet divergence, which can be properly treated by 
renormalization. This must be done in such a way as 
to maintain the correct underlying two-body resonance 
physics for any momentum cutoff K in the field the- 
ory. One should consider the delta function interactions 




FIG. 1. Binding energies as a function of magnetic field. 
The solid line is the coupled channels result for the most 
accurate rubidium interactions. This is compared with the 
binding energy resulting from a contact scattering model with 
<2b g = — 450ao (where ao is the Bohr radius), A^i = 2.23/is 
(with fin the Bohr magneton), and go = 3.11 x 10~ 38 Jm 3 ' 2 
(open squares). 

The outcome of the experiment of Donley et al. S 
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closely resembles the seminal experiments on Ramsey 
fringes in atomic beam physics [jl3j. The starting point 
is a condensate of 85 Rb atoms in the = 2,— 2} 

state, at a magnetic field where the scattering length is 
close to zero. Two magnetic field pulses are applied (a 
schematic is shown in Fig. |2|), each of which brings the 
condensate close to resonance. The two pulses are sep- 
arated by a free evolution interval t cvo i vc , during which 
time the magnetic field is increased to move the system 
further away from resonance. After this pulse sequence, 
the remaining number of atoms in the condensate is mea- 
sured, which is then called the remnant. Also a burst of 
noncondensate atoms is observed. The populations of the 
remnant and the burst both show oscillations as function 
of tcvoivc at a frequency that corresponds to the binding 
energy of the molecular state at the intermediate field. 
The sum of the populations of the remnant and burst 
do not add up to the initial number, implying a miss- 
ing component. The theoretical solutions which follow 
in this Letter can be directly compared with the results 
presented in the figures of Donley et al. [§. 
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FIG. 2. A typical magnetic field pulse sequence as a func- 
tion of time. The interval i ovo i ve is modified in the experiment. 
The position of the resonance is indicated by the dashed line. 

In spite of the fact that the mean-field Eqns. (||)-(0) 
were derived for a homogeneous system, we may apply 
our theory to a trapped gas. At the highest energy scales 
reached during the evolution, the velocity of the atoms 
is sufficient to move only one hundredth of the oscillator 
ground state size in the tightest direction during the full 
evolution time. Therefore, a complete quantum descrip- 
tion of the oscillator levels is not required. Instead, we 
use a local density approximation and perform a Gaus- 
sian average over the densities of the gas. For each value 
of the density, we solve the time-dependent equations 
modifying the detuning according to the time-dependent 
field given in Fig. ||. In Fig. || we show a time evolution 
of the atomic condensate, and the corresponding time 
evolution of the molecular condensate. At the end of the 



first pulse about 25% of the condensate atoms have been 
converted into other components. The growth of pop- 
ulation in the molecular condensate takes place mostly 
during the final ramp. It is notable that the fraction of 
molecules does not account for the missing atoms. In 
fact, the atoms are mainly transfered to the normal and 
anomalous densities which are ascribed to the noncon- 
densate component. 
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FIG. 3. Fraction of the atomic condensate as a function 
of time. We identify these atoms as the remnant of Donley et 
al. The calculation is done for the pulse sequence given 
in Fig. |2|, for a density n = 3.9 x 10 12 cm~ 3 . The inset shows 
the analogous graph for the molecular condensate. 

The growth of the noncondensate component can been 
seen in Fig. [|. The function Gat(Y = 0) represents the 
density of noncondensate atoms, which can be seen to 
oscillate out of phase with the atomic condensate dur- 
ing fevoive- It is remarkable that in spite of the fact that 
the Ramsey fringes occur at the molecular binding energy 
with significant visibility, the population of the molecular 
condensate remains small. A much larger fraction is con- 
verted into strongly correlated atom pairs, encapsulated 
by the normal and anomalous densities. Interestingly, the 
anomalous density is the same aspect of the field theory 
which accounts for Cooper pairing in a nonidcal Fermi 
gas and gives rise to superfluidity at temperatures be- 
low the critical value. The explanation for the observed 
growth of the pairing field rather than the molecular con- 
densate is due to the close proximity of the bound state 
to threshold. The range of the molecular bound state 
stretches in this case to very large internuclear distances, 
something which has much more overlap with the delo- 
calized Ga pairing field than with the localized closed 
channel state. 

In Fig. |B| we show the population of atomic conden- 
sate and noncondensate atoms, obtained at the end of 
the pulse sequence, as a function of the evolution time 
^evolve- The sum of these two numbers (squares) equals 
the total number of initial condensate atoms minus twice 
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the number of molecular condensate atoms due to parti- 
cle conservation. The frequency of the oscillations agrees 
with the binding energy of the renormalized potential, 
given by the open squares in Fig. [j]. When we com- 
pare these curves with Fig. 5 of Donley et at M, we see 
that we can clearly identify the remnant observed in the 
experiment as the atomic condensate component of the 
quantum field theory. The experimental data closely re- 
sembles the solid curve both in offset and in amplitude. 
Similarly, the noncondensate atoms can be identified as 
the burst atoms. Since Ga/-(0 is a correlation function, 
it is straightforward to determine the energy of the non- 
condensed atoms which are produced. This manifests 
as the spatial decay rate of the correlation function in 
the r-direction, which can be converted into an average 
energy. The energy which results is comparable to the 
experimentally determined energy range for the burst. 
The missing atoms are also elucidated. Note that there 
is a large time interval between the final ramp and the 
time at which the atomic absorption imaging takes place. 
These weakly bound molecules may decay to lower vibra- 
tional states via a collision with a third atom, resulting in 
large kinetic energies for both scattering partners. Such 
atoms would not be observed. 




FIG. 4. Normal field Gjv(j") as a function of time and 
distance for the same calculation illustrated in Fig. ^. We 
identify these atoms as the burst atoms of Donley et al. Q. 
The spread in r-space is a measure of the energy of these 
noncondensate atoms. 

We have repeated our calculation for a different exper- 
imental situation Jl4| with a factor 10 larger density, and 
a different time-dependence of the field. Here the number 
of remnant atoms is lower than the burst atoms, so that 
the position of the fringes shown in Fig. [| are switched 
and show indications of damping. We again get good 
agreement with the experiment, and the appearance of 
damping of the Ramsey fringes in the theory is due to 
density-dependent inhomogeneous dephasing. This ef- 
fect is due to a relatively small shift in the oscillation 



frequency that shows a strong dependence on density. 
Finally we note that both theory and experimental data 
exhibit a notable phase-shift between the position of the 
fringes of the atomic condensate and the noncondensate 
atoms. 
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FIG. 5. Oscillations between the atomic condensate (solid 
line) and the normal field Gn(0) (dashed line). These two 
numbers add up to the total number of recovered atoms 
(squares), which excludes the molecular component. The cal- 
culation is performed for a mean density of (n) = 3.9 x 10 12 
cm -3 . This case can be directly compared with Fig. 6 of 
Donley et al. 

In conclusion, we have used a resonance effective field 
theory which includes an accurate description of the two- 
body bound state and scattering physics to describe a 
recent experiment at JILA. Because of the complexity 
of the mean-field physics relevant to Feshbach resonance 
scattering, it has not been possible previously to pro- 
vide this kind of quantitative comparison. We are able 
to unambiguously identify the observed remnant and 
burst. The pairing field associated with the nonconden- 
sate atoms plays a crucial role in our calculated evolu- 
tion. This pairing field is analogous to the formation of 
Cooper pairs in a superfluid Fermi gas. The ability to 
determine the coupling constants from known two-body 
rubidium physics allows us to make these comparisons 
with experimental data with no adjustable parameters. 
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